Design and Implementation of a High-Speed Readout and Control System for a Digital Tracking Calorimeter for proton CT

Particle therapy, a non-invasive technique for treating cancer using protons and light ions, has become more and more common. For example, a particle treatment facility is currently being built, in Bergen, Norway. Proton beams deposit a large fraction of their energy at the end of their paths, i.e., the delivered dose can be focused on the tumor, sparing nearby tissue with a low entry and almost no exit dose. A novel imaging modality using protons promises to overcome some limitations of particle therapy and allowing to fully exploit its potential. Being able to position the so-called Bragg peak accurately inside the tumor is a major advantage of charged particles, but incomplete knowledge about a crucial tissue property, the stopping power, limits its precision. A proton CT scanner provides direct information about the stopping power. It has the potential to reduce range uncertainties significantly, but no proton CT system has yet been shown to be suitable for clinical use. The aim of the Bergen proton CT project is to design and build a proton CT scanner that overcomes most of the critical limitations of the currently existing prototypes and which can be operated in clinical settings. A proton CT prototype, the Digital Tracking Calorimeter, is being developed as a range telescope consisting of high-granularity pixel sensors. The prototype is a combined position-sensitive detector and residual energy-range detector which will allow a substantial rate of protons, speeding up the imaging process. The detector is single-sided, meaning that it employs information from the beam delivery system to omit tracker layers in front of the phantom. The detector operates by tracking the charged particles traversing through the detector material behind the phantom. The proton CT prototype will be used to determine the feasibility of using proton CT to increase the dose planning accuracy for particle treatment of cancer cells. The detector is designed as a telescope of 43 layers of sensors, where the two front layers act as the position-sensitive detector providing an accurate vector of each incoming particle. The remaining layers are used to measure the residual energy of each particle by observing in which layer they stop and by using the cluster size in each layer. The Digital Tracking Calorimeter employs the ALPIDE sensor, a monolithic active pixel sensor, each utilizing a 1.2Gb/s data link. Each layer of 18×27 cm consists of 108 ALPIDE sensors, roughly corresponding to the width and height of the head of a grown person. The sensors are connected to intermediary transition boards that route the data and control links to dedicated readout electronics and supply the sensors with power. The readout unit is the main component of both the data acquisition and the detector control system. The power control unit controls the power supply and monitors the current usage of the sensors. Both of these devices are mainly implemented in FPGAs. The main purpose of this work has been to explore and implement possible design solutions for the proton CT electronics, including the front-end, as well as the readout electronics architecture. The resulting architecture is modular, allowing the further scale-up of the system in the future. A major obstacle to the design is the high amount of sensors and the corresponding high-speed data links. Thus, a large emphasis has been on the signal integrity of the front-end electronics and a dynamic phase alignment sampling method of the readout electronics firmware. The readout FPGA employs regular I/O pins for the high-speed data interface, instead of high-speed transceiver pins, which significantly reduces the magnitude of the data acquisition system. A consistent design approach with detailed and systematic verification of the FPGA firmware modules, along with a continuous integration build system, has resulted in a stable and highly adaptive system. Significant effort has been put into the testing of the various system components. This also includes the design and implementation of a set of production test tools for use during the manufacturing of the detector front-end.Doktorgradsavhandlin

Image quality of list-mode proton imaging without front trackers

List mode proton imaging relies on accurate reconstruction of the proton most likely path (MLP) through the patient. This typically requires two sets of position sensitive detector systems, one upstream (front) and one downstream (rear) of the patient. However, for a clinical implementation it can be preferable to omit the front trackers (single-sided proton imaging). For such a system, the MLP can be computed from information available through the beam delivery system and the remaining rear tracker set. In this work, we use Monte Carlo simulations to compare a conventional double-sided (using both front and rear detector systems) with a single-sided system (only rear detector system) by evaluating the spatial resolution of proton radiographs (pRad) and proton CT images (pCT) acquired with these set-ups. Both the pencil beam spot size, as well as the spacing between spots was also adjusted to identify the impact of these beam parameters on the image quality. Relying only on the pencil beam central position for computing the MLP resulted in severe image artifacts both in pRad and pCT. Using the recently extended-MLP formalism that incorporate pencil beam uncertainty removed these image artifacts. However, using a more focused pencil beam with this algorithm induced image artifacts when the spot spacing was the same as the beam spot size. The spatial resolution tested with a sharp edge gradient technique was reduced by 40% for single-sided (MTF10% = 3.0 lp/cm) compared to double-sided (MTF10% = 4.9 lp/cm) pRad with ideal tracking detectors. Using realistic trackers the difference decreased to 30%, with MTF10% of 4.0 lp/cm for the realistic double-sided and 2.7 lp/cm for the realistic single-sided setup. When studying an anthropomorphic paediatric head phantom both single- and double-sided set-ups performed similarly where the difference in water equivalent thickness (WET) between the two set-ups were less than 0.01 mm in homogeneous areas of the head. Larger discrepancies between the two set-ups were visible in high density gradients like the facial structures. A complete CT reconstruction of a Catphan$^{\circledR}$ module was performed. Assuming ideal detectors, the obtained spatial resolution was 5.1 lp/cm for double-sided and 3.8 lp/cm for the single-sided setup. Double- and single-sided pRad with realistic tracker properties returned a spatial resolution of 3.8 lp/cm and 3.2 lp/cm, respectively. Future studies should investigate the development of dedicated reconstruction algorithms targeted for single-sided particle imaging.publishedVersio

Production of charged pions, kaons, and (anti-)protons in Pb-Pb and inelastic pp collisions at √sNN=5.02 TeV

Midrapidity production of π±, K±, and (¯p)p measured by the ALICE experiment at the CERN Large Hadron Collider, in Pb-Pb and inelastic pp collisions at √sNN=5.02 TeV, is presented. The invariant yields are measured over a wide transverse momentum (pT) range from hundreds of MeV/c up to 20 GeV/c. The results in Pb-Pb collisions are presented as a function of the collision centrality, in the range 0–90%. The comparison of the pT-integrated particle ratios, i.e., proton-to-pion (p/π) and kaon-to-pion (K/π) ratios, with similar measurements in Pb-Pb collisions at √sNN=2.76 TeV show no significant energy dependence. Blast-wave fits of the pT spectra indicate that in the most central collisions radial flow is slightly larger at 5.02 TeV with respect to 2.76 TeV. Particle ratios (p/π, K/π) as a function of pT show pronounced maxima at pT≈3 GeV/c in central Pb-Pb collisions. At high pT, particle ratios at 5.02 TeV are similar to those measured in pp collisions at the same energy and in Pb-Pb collisions at √sNN=2.76 TeV. Using the pp reference spectra measured at the same collision energy of 5.02 TeV, the nuclear modification factors for the different particle species are derived. Within uncertainties, the nuclear modification factor is particle species independent for high pT and compatible with measurements at √sNN=2.76 TeV. The results are compared to state-of-the-art model calculations, which are found to describe the observed trends satisfactorily.publishedVersio

J/ψ elliptic and triangular flow in Pb-Pb collisions at √sNN = 5.02 TeV

The inclusive J/ψ elliptic (v2) and triangular (v3) flow coefficients measured at forward rapidity (2.5 < y < 4) and the v2 measured at midrapidity (|y| < 0.9) in Pb-Pb collisions at sNN−−−√ = 5.02 TeV using the ALICE detector at the LHC are reported. The entire Pb-Pb data sample collected during Run 2 is employed, amounting to an integrated luminosity of 750 μb−1 at forward rapidity and 93 μb−1 at midrapidity. The results are obtained using the scalar product method and are reported as a function of transverse momentum pT and collision centrality. At midrapidity, the J/ψ v2 is in agreement with the forward rapidity measurement. The centrality averaged results indicate a positive J/ψ v3 with a significance of more than 5σ at forward rapidity in the pT range 2 < pT < 5 GeV/c. The forward rapidity v2, v3, and v3/v2 results at low and intermediate pT (pT ≲ 8 GeV/c) exhibit a mass hierarchy when compared to pions and D mesons, while converging into a species-independent curve at higher pT. At low and intermediate pT, the results could be interpreted in terms of a later thermalization of charm quarks compared to light quarks, while at high pT, path-length dependent effects seem to dominate. The J/ψ v2 measurements are further compared to a microscopic transport model calculation. Using a simplified extension of the quark scaling approach involving both light and charm quark flow components, it is shown that the D-meson vn measurements can be described based on those for charged pions and J/ψ flow.publishedVersio

Measurement of electrons from heavy-flavour hadron decays as a function of multiplicity in p-Pb collisions at √sNN = 5.02 TeV

The multiplicity dependence of electron production from heavy-flavour hadron decays as a function of transverse momentum was measured in p-Pb collisions at √sNN = 5.02 TeV using the ALICE detector at the LHC. The measurement was performed in the centre-of-mass rapidity interval −1.07 < ycms< 0.14 and transverse momentum interval 2 < pT< 16 GeV/c. The multiplicity dependence of the production of electrons from heavy-flavour hadron decays was studied by comparing the pT spectra measured for different multiplicity classes with those measured in pp collisions (QpPb) and in peripheral p-Pb collisions (Qcp). The QpPb results obtained are consistent with unity within uncertainties in the measured pT interval and event classes. This indicates that heavy-flavour decay electron production is consistent with binary scaling and independent of the geometry of the collision system. Additionally, the results suggest that cold nuclear matter effects are negligible within uncertainties, in the production of heavy-flavour decay electrons at midrapidity in p-Pb collisions.publishedVersio

Measurement of strange baryon–antibaryon interactions with femtoscopic correlations

Two-particle correlation functions were measured for pp, p, p, and pairs in Pb–Pb collisions at √sNN = 2.76 TeV and √sNN = 5.02 TeV recorded by the ALICE detector. From a simultaneous fit to all obtained correlation functions, real and imaginary components of the scattering lengths, as well as the effective ranges, were extracted for combined p and p pairs and, for the first time, for pairs. Effective averaged scattering parameters for heavier baryon–antibaryon pairs, not measured directly, are also provided. The results reveal similarly strong interaction between measured baryon–antibaryon pairs, suggesting that they all annihilate in the same manner at the same pair relative momentum k∗. Moreover, the reported significant non-zero imaginary part and negative real part of the scattering length provide motivation for future baryon–antibaryon bound state searches.publishedVersio

Measurement of isolated photon-hadron correlations in √sNN = 5.02 TeV pp and p-Pb collisions

This paper presents isolated photon-hadron correlations using pp and p-Pb data collected by the ALICE detector at the LHC. For photons with |η|<0.67 and 12<pT<40 GeV/c, the associated yield of charged particles in the range |η|<0.80 and 0.5<pT<10 GeV/c is presented. These momenta are much lower than previous measurements at the LHC. No significant difference between pp and p-Pb is observed, with Pythia 8.2 describing both data sets within uncertainties. This measurement constrains nuclear effects on the parton fragmentation in p-Pb collisions, and provides a benchmark for future studies of Pb-Pb collisions.publishedVersio

Constraining the Chiral Magnetic Effect with charge-dependent azimuthal correlations in Pb-Pb collisions at sNN\sqrt{s_{NN}} = 2.76 and 5.02 TeV

Systematic studies of charge-dependent two- and three-particle correlations in Pb-Pb collisions at $\sqrt{s_{NN}}$ = 2.76 and 5.02 TeV used to probe the Chiral Magnetic Effect (CME) are presented. These measurements are performed for charged particles in the pseudorapidity (η) and transverse momentum (pT) ranges |η| < 0.8 and 0.2 < pT < 5 GeV/c. A significant charge-dependent signal that becomes more pronounced for peripheral collisions is reported for the CME-sensitive correlators γ1, 1 = 〈cos(φα + φβ − 2Ψ2)〉 and γ1, − 3 = 〈cos(φα − 3φβ + 2Ψ2)〉. The results are used to estimate the contribution of background effects, associated with local charge conservation coupled to anisotropic flow modulations, to measurements of the CME. A blast-wave parametrisation that incorporates local charge conservation tuned to reproduce the centrality dependent background effects is not able to fully describe the measured γ1,1. Finally, the charge and centrality dependence of mixed-harmonics three-particle correlations, of the form γ1, 2 = 〈cos(φα + 2φβ − 3Ψ3)〉, which are insensitive to the CME signal, verify again that background contributions dominate the measurement of γ1,1.publishedVersio

Multiplicity dependence of light (anti-)nuclei production in p–Pb collisions at √sNN=5.02 TeV

The measurement of the deuteron and anti-deuteron production in the rapidity range −1 < y < 0 as a function of transverse momentum and event multiplicity in p–Pb collisions at √sNN = 5.02 TeV is presented. (Anti-)deuterons are identified via their specific energy loss dE/dx and via their time-offlight. Their production in p–Pb collisions is compared to pp and Pb–Pb collisions and is discussed within the context of thermal and coalescence models. The ratio of integrated yields of deuterons to protons (d/p) shows a significant increase as a function of the charged-particle multiplicity of the event starting from values similar to those observed in pp collisions at low multiplicities and approaching those observed in Pb–Pb collisions at high multiplicities. The mean transverse particle momenta are extracted from the deuteron spectra and the values are similar to those obtained for p and particles. Thus, deuteron spectra do not follow mass ordering. This behaviour is in contrast to the trend observed for non-composite particles in p–Pb collisions. In addition, the production of the rare 3He and 3He nuclei has been studied. The spectrum corresponding to all non-single diffractive p-Pb collisions is obtained in the rapidity window −1 < y < 0 and the pT-integrated yield dN/dy is extracted. It is found that the yields of protons, deuterons, and 3He, normalised by the spin degeneracy factor, follow an exponential decrease with mass number.publishedVersio

J/ψ production as a function of charged-particle multiplicity in p-Pb collisions at √sNN = 8.16 TeV

Inclusive J/ψ yields and average transverse momenta in p-Pb collisions at a center-of-mass energy per nucleon pair sNN−−−√ = 8.16 TeV are measured as a function of the charged-particle pseudorapidity density with ALICE. The J/ψ mesons are reconstructed at forward (2.03 < ycms < 3.53) and backward (−4.46 < ycms < −2.96) center-of-mass rapidity in their dimuon decay channel while the charged-particle pseudorapidity density is measured around midrapidity. The J/ψ yields at forward and backward rapidity normalized to their respective average values increase with the normalized charged-particle pseudorapidity density, the former showing a weaker increase than the latter. The normalized average transverse momenta at forward and backward rapidity manifest a steady increase from low to high charged-particle pseudorapidity density with a saturation beyond the average value.publishedVersio